Anderson transition in one-dimensional systems with spatial 

disorder 

Rabah Benhenni 1 , Khaled Senouci 1,2 J*| Rachid Bouamrane 1 , and Nouredine Zekri 1 
1 U.S.T.O., Departement de Physique L.E.P.M., 
B. P. 1505 El M'Naouar, Oran, Algeria and 
2 Universite de Mostaganem Abdelhamid Ibn Badis, 
Departement de Physique, B.P. 227, 
Route Belhacel, 27000, Mostaganem, Algeria 

Abstract 

A simple Kronig-Penney model for one-dimensional (ID) mesoscopic systems with 5 peak poten- 
tials is used to study numerically the influence of a spatial disorder on the conductance fluctuations 
and distribution at different regimes. We use the Levy laws to investigate the statistical properties 
of the eigenstates. We found the possibility of an Anderson transition even in ID meaning that the 
disorder can also provide constructive quantum interferences. We found at this transition that the 
conductance probability distribution has a system-size independent shape with large fluctuations in 
good agreement with previous works. In these ID systems, the metallic phase is well characterized 
by a Gaussian conductance distribution. Indeed, the results for the conductance distribution are 
in good agreement with the previous works in 2D and 3D systems for other models. 
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I. INTRODUCTION 



The disorder-induced metal-insulator transition (MIT) has been studied extensively for 

n n n 

decades |2|, [3J and continues to attract much attention. Scaling theory |2J predicts that 
all eigenstates of noninteracting electrons are localized in one-dimensional (ID) and two- 
dimensional (2D) systems for any amount of disorder. It is commonly believed that MIT 
occurs only for dimensions d > 2 and the system remain insulator for d < 2 [4|. However, 
recently, it was suggested that a 2D Anderson model of localization with purely off-diagonal 
disorder might violate this general statement since non-localized states were found at the 
band center [5|. It was found that a transition present in this 2D model can be described by 
a localization length which diverges with a power-law behaviour. The possible existence of 
the Anderson transition in 2D systems without interaction and spin-orbit effects becomes 
recently a subject of controversy in the literature 

a a. 

More recently, Asada et al. 
studied the f3 function that describes the scaling of the quantity A as: 



P(Ln A) 



dLn A 



(1) 



dLn L 

where A is the ratio of the quasi-one-dimensional localization length to the system width 
for electrons on a long quasi-one-dimensional system of width L. They indicated the possibil- 
ity of an Anderson transition for dimensions d < 2 in disordered systems of non-interacting 
electrons. On the other hand, strong numerical evidence of a mobility edge was found in 
disordered photonic systems in two dimensions l| . Recently, surp rising results were found 
in the properties of disordered graphene systems [h], 11, 12, 131]. The latest studies, both 



theoretical and experimental, led to the amazing conclusion that there is no localization in 
disordered graphene, even in the one-dimensional (ID) situation. 



It is well known that the conductance g is not a self-averaged quantity [14[ and therefore 
fluctuates with the Fermi energy, chemical potential and sample size. Since the conductance 
does not obey the central limit theorem [15|. it is necessary to investigate not only the 
first two moments but the entire probability distribution. Numerical results in 2D and 
3D disordered systems showed that the conductance is normally distributed in the metallic 



regime whi 
was found 



e for strongly localized systems (insulating regime) a log-normal distribution 



161 ]. The exact form of the probability distribution at the transition is not well 
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known. In such a regime, it was proven that the conductance distribution is independent 
of the microscopic details of the model (determined by the distribution of the disorder), 
of the system size and of the position of the critical point which separates the metallic 
and the localized regime in the phase space of external parameters (energy, disorder). The 
universality of the conductance distribution was studied and confirmed by Markos et al. for 

18) . The system-size invariance of the conductance distribution at 



2D and 3D models 
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the critical points of the MIT was confirmed for 3D and 4D systems in 

The aim of this paper is to investigate the possibility of an Anderson transition in ID 
spatially disordered systems of noninteracting electron. Three regimes of electron transport 
are studied: the insulating regime corresponding to strong disorder, the metallic regime 
(corresponding to an infinitesimal disorder) and the transition regime. In the present work, 
Levy statistics B] are used. 



II. MODEL DESCRIPTION 

We consider a Kronig-Penney model applied to a ID system of delta potentials with 
random positions (spatial disorder) . The corresponding Schrodinger equation can be read: 



d 2 



dx 2 



+ I 3 S ( X ~ X n) \ ®( X ) = E ^( X ) 



(2) 



Here \l/(x) is the single particle wave- function at x, (3 the potential strength and E the 
single particle energy in units of h 2 /2m with m the electronic effective mass. The two ends 
of the system are assumed to be connected ohmically to ideal leads (where the electron 
moves freely) 

The second order differential equation (2) can be mapped by means the Poincare map 
representation 23]: 



cos(fc/ n+ i) + 



sin(/c/ n 



cos (£;/„) + f3 n 



sin(/c/ n+ i) 
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where \l/ n is the value of the wave-function at site n and k 
number l n+ i = x n +i — x n is the inter-atomic spacing. l n+ i 



n+lj 



n-1 



(3) 



sm(kl n ) 

E is the electron wave 
a + s, s being a random 



variable uniformly distributed as — W/2 < s < W/2 (W being the degree of the disorder). 



For equally spaced potentials, the interatomic spacing l n+ \ 
equation (2) reduces to that found by Soukoulis et al. [2i 



l n = a (W = 0) and 



n-l 



2 cos VE + p n 



sin v i£ 



(4) 



The solution of equation (3) is carried out iteratively by taking the two initial wave 
functions at sites 1 and 2 : = exp(— ik) and ^2 = exp(— 2ik). We consider here an 
electron having a wave number k incident at site iV + 3 from the right (by taking the chain 
length L = N, i.e. N + 1 scatterers). The transmission coefficient (T) reads 



T 



exp(—2ikl n ) 



(5) 



- * N +3exp(-ikl n )\ 2 

The dimensionless conductance (g = ) can be obtained from the transmission coeffi 
cient T via the Landauer formula for ID systems [24] : 



IT 
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1 — T 

where the factor two arises from the two possible states of the electron spin, 
and the variance of conductance (Ln(g)) reads: 



(6) 



var(Ln(g)) = < Ln(g) 2 > - < Ln(g) > 2 (7) 
where < ... > denotes an average over different realizations of the disorder. 



III. RESULTS AND DISCUSSION 



In this section, we discuss numerical results of the conductance probability distribution 
for different amount of spatial disorder of ID mesoscopic systems of delta potentials. In order 
to obtain the probability distribution of the conductance, we build a statistical ensemble of 
10 4 samples differing only in the realization of the disorder. 

Firstly, we examine the effect of this kind of disorder on electronic eigenstates by calcu- 
lating the transmission coefficient and the conductance probability distribution. Figure 1 
shows the scaling of the transmission coefficient < —Ln(T) > for different degrees of spatial 
disorder and for two values of electron energy E = 4 (corresponding to an energy in the 
band gap of the periodic system) and E = 12 (corresponding to an energy in the allowed 
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E = 4 


E = 4 


E = 12 


E = 12 
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c 


7 


0.005 


4797.866 


2.084 x 10~ 4 


1.7832 


0.5607 


0.1 


1205.313 


8.2966 x 10~ 4 


1.8082 


0.5530 


0.2 


0301.204 


0.00332 


1.9205 


0.52068 



TABLE I: Calculated results of the Lyapunov exponent 7 and the localization length £ as a function 
of the disorder W. 



band, see inset of Figure l.a). For strong disorder the transmission coefficient decreases ex- 
ponentially with the length scale. Furthermore, the localization length ( and the Lyapunov 
exponent 7 = l/£(slope of the curve) are deduced. The calculated results are presented in 
table I. Fig. la shows that increasing disorder leads to strong localization (since the local- 
ization length ( decreases with the disorder strength W). This is in good agreement with 
previous results with topological disorder i.e the disorder localizes the electronic eigenstates. 
However, this is not the case as shown in Fig. lb where ( increases with the disorder. This 
behaviour was previously pointed out by Nimour et al. for ID spatially disordered systems 
with finite width potentials 25[. More recently, it has been shown that the disorder could 
either suppress or enhance the transmission in disordered grapheme superlattices 26]. 



Since our system is disordered, its conductance depends on specific realizations of dis- 
order (for a given L), it is then appropriate to study the whole probability distribution of 
conductance. For strong disorder where w = 0.2, we have ( = 300.2, which is smaller than 
the size of the system. The localization of the electronic eigenstates is confirmed in figure 
2a where the probability distribution of the natural logarithm of the conductance (Ln(g)) 
is plotted. In this figure the probality distribution is Gaussian indicating the general be- 
haviour in this regime This figure shows also that the mean conductance decreases 



as the system size increases. In a previous interesting paper, Vagner et al. 27] studied 
analytically electron transport in a one-dimensional wire with disorder modeled as a chain 
of randomly positioned scatterers with 5-shaped impurity potential (similar to our model). 
They found that the distribution P(f) of the variable / = ln(l + p) (p being the resistance) 
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has a non-Gaussian behaviour in the limit of weak disorder. This result is confirmed for 
the conductance in our system. It is clearly seen from figure 2b that when the disorder de- 
creases, the conductance distribution shows a deviation from its log-normal form indicating 
derealization of the eigenstates. 

In order to further understand the behaviour of the conductance, we investigate the 
conductance distribution P(g) for different amounts of spatial disorder and for different 
system sizes. All these distributions show long power law tails decreasing for large values 
of the conductance g (see Figure 3). Therefore, we use the Levy statistics L^(Z) of index 
H which decreases as Z~( 1+ ^ for large values of Z 22|. It is found that the conductance 
distribution P(g) behaves as g~( 1+ ^ for large values of g. The exponent \i is then extracted 
from the log-log plot of P(g) for large values of g which is linear with a slope equal to — 
(see the inset). It is known that if fi > 2, the probability distribution is normal. On the 
other hand, if /i < 2 it means that the distribution is log-normal. In figure 4, the index // is 
plotted as a function of degree of spatial disorder W and for different system size L. In this 
figure two distinguishable regions are shwn: a region of fi < 2 (corresponding to an insulating 
regime) and a region of fi > 2 (corresponding to a metallic regime for small disorder W). 
The intersection point of curve fi(W) with the strait line /i = 2 corresponds to the metal- 
insulator transition region. The critical point separating the metallic and insulating regimes 
corresponds to a critical disorder with W c = 0.0035. Figure 4 shows also that the index /i is 
independent of the system size at the transition. 

Let us now examine the conductance probability distribution for different amounts of 
spatial disorder. In Figure 5 are plotted the probability distributions of Ln(g) (Fig.5a) 
and g (Fig.5b) for different system size (L = 500, 700, 800 and 900) for the critical disorder 
W c = 0.0035. The conductance distribution seems to be neither normal nor log-normal. 
The size independence of the conductance distribution is in agreement with the observed 

nnn 

one for 2D and 3D systems [13, H9J, |21[ at the transition. The long tail of the distribution in 
Tig. 5a is representative for large fluctuations in good agreement with the results of Shapiro 
28| and Shapiro and Cohen 29] for the metal-insulator transition. 

For an infinitesimal disorder W < W c corresponding to fi > 2, where w = 0.0025, we 
have ( = 47234, which is much larger than the system size L. The conductance distribution 
has a Gaussian form indicating the metallic regime (Figure 6). This figure shows the general 
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behaviour for this regime. Once the system is away from the critical point, the conductance 
distribution P(g) begins to show size dependence. Distribution P(g) moves toward higher 
values of the conductance as the system size increases. 



IV. CONCLUSION 



We have used the Kronig-Penney model in a simple one-dimensional (ID) system with 
spatial disorder to determine the metal-insulator transition and examine the size indepen- 
dence of its distribution at this transition. To find the critical point, we have used stable 
Levy laws 22j. We found the possibility of an Anderson transition even in ID meaning 
that the disorder can also provide constructive quantum interferences. We found the critical 
disorder W c for this transition. Indeed, the results are in good agreement with the previous 



works in 2D and 3D systems for other models 13, ll9|, |21| for metal-insulator transition. The 



conductance probability distribution is found to be Gaussian in the metallic regime. At the 
transition, the distribution P(lng) exhibits the typical asymmetric behaviour while in the 
localized regime, the distribution is log normal. It is important to study the universality of 
the conductance distribution at the transition in a system of finite width potentials where 
the conductance fluctuations are less important in comparison to the present model 30], 
and to find the critical points of the transition in the energy- disorder phase space and for 
different kinds of disorder. It is also important to investigate the effect of spatial disorder 
on conductance fluctuations in different regimes. These problems will be the subject of a 
forthcoming paper. 
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Figure captions 

Figure 1: Transmission coefficient <-Ln(T)> as a function of system size averaged over 
1000 realizations of the same system for V = 2 and different amounts of spatial disorder W 
and for a) E — 4 (in the miniband). b) E — 12 (in the gap). Inset in a: Transmission Vs 
energy for a periodic system. 

Figure 2: Probability distribution of —log(g) for E = 4, V = 2, L = 1600 and for a) 
W = 0,2. b) W = 0,0036. 

Figure 3: Conductance probability distribution for L = 1500, E = 4 and W = 0.0034. 
Inset: log-log plot of the tail of the distribution, the strait line corresponds to the best 
power-law fit to the tail P(g) ~ g~( 1+ v) with fi = 1.02. 

Figure 4: Levy exponent \i as a function of disorder for E = 4, V = 2 and different 
system size L. 

Figure 5: a) Probability distribution of the conductance g for E = 4 and V = 2. 
b) Probability distribution of — ln(g) compared with a Gaussian with the same mean and 
variance. 

Figure 6: Conductance probability distribution in the metallic regime for different sys- 
tem size L with W = 0.0025, E = 4 and V = 2 compared with a Gaussian with the same 
mean and variance. 
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